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Abstract 

We examine the sensitivity of wavefunction intensities in chaotic quantum 
systems to small changes in an arbitrary external perturbation. A universal 
scaling is proposed for all three Dyson ensembles and a novel theoretical 
approach is used to determine exact expressions for systems which violate 
T-invariance. Analytical results are compared with numerical simulations of 
tight-binding Anderson Hamiltonians. 
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Systems which are classically chaotic show a wide degree of universality in their quantum 
properties The generality of this description unites a large class of systems which 

include examples as diverse as the spectra of atomic nuclei ||, Rydberg atoms in strong 
magnetic fields |4[], weakly disordered metallic grains |5|J|, and ballistic electronic cavities 
0. Studies have focussed largely on the statistical properties of spectra and wavefunction 
intensities. Recently, there has been an attempt to extend the universal description to 
encompass the response of spectra to changes in an arbitrary external perturbation |S],|S]]. In 
this letter we will show how this description extends to the study of parametric correlations 
in properties of wavefunctions amplitudes. 

Traditional studies have involved defining spectral averages at fixed values of external 
fields ||10|| . From an experimental point of view it is clearly more desirable to understand 
what happens to the statistics of a given eigenfunction at a certain position upon varying an 
external parameter such as magnetic field. Such information is contained within the joint 
distribution function 

W(v,v;X) = (5(v + v/2-n\il> v (r,X + X/2)| 2 ) 

x6(v - v/2 - fi|W(r, X - X/2)| 2 )), (1) 

where ip v {r,X) represents the eigenfunction of the v-t\i state at a position r, say inside 
a quantum dot, X parametrizes some external perturbation such as magnetic field, and 
fl denotes the total volume of the system. Statistical averaging, denoted by (• • •) can in 
principle be performed over some range of energy, parameter, or over the statistical ensemble 
if one is defined. Our study is motivated by a number of recent experiments that, in principle, 
allow W to be measured directly. These include STM measurements in nanostructures or 
"quantum corrals" |TT], microwave cavities [I2|, and small acoustic resonators [[13 . 



The universality of the statistical properties of weakly disordered metallic grains has been 
verified through the supersymmetry approach [0]. On time scales longer than the typical 
diffusion time = L 2 /D, where D = v\rjd denotes the classical diffusion coefficient, 
statistical properties of the spectra and wavefunctions coincide with those of random matrix 
ensembles. In this regime there exist only two relevant parameters, the average level spacing 
A = (E u+ i — E u ), and the mean square level gradient C(0) = ((de u /dX) 2 ) (typically 
C(0) ~ l/tflA). The dependence of statistical properties on both can be removed by a 
straightforward unfolding of the energy levels e M = E^/A, and a rescaling of the perturbation 



x = XyC(O). Strong numerical evidence (for a review, see Refs. H|T4|]) suggests that this 
type of universality extends to all non-integrable systems and serves as a characterization 
of the general phenomenum of quantum chaos. 

Below we will show that, asymptotically in the limit x — > 0, the joint distribution function 
takes the universal form 

W(v, v; x) = P p (v) - — -^e B p (-^7=) , (2) 

where P p (v) = (vf3 /2) /3 / 2 - 1 e~ s/3/2 (/3/2)/r(/3/2) denotes the \ 2 distribution of wavefunction 
intensities for Dyson ensembles of Hermitian random matrices, and Bp(s) is a universal 
function whose form depends only on the symmetries of the system. By convention, T- 
invariant ensembles that conserve/ violate spin rotational symmetry are termed orthogonal 
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(P = l)/symplectic (f3 = 4) respectively. Those which violate T-invariance are known as 
unitary (J3 = 2). For the latter, a new approach to calculate the statistical properties of 
spectral determinants [[15] yields the following expression for the correlator 

35 + 14s 2 + 3s* 

The powerlaw behavior of s ^> 1, which implies a large sensitivity, is the related to the 
algebraic level repulsion of the Wigner-Dyson statistics. 

Before comparing W with numerical simulations, we will show how random matrix theory 
(RMT) leads to the general scaling in Eq. (fj). We emphasize that the same conclusions could 
be drawn from the study of a microscopic theory involving weakly disordered metallic grains. 
Consider the Hamiltonian 

H = H + XH 1 , (4) 

where Hq and Hi denote N x N complex random Hermitian matrices which belong to the 
Gaussian unitary ensemble defined by the distribution P{H) dH oc exp[— (N/2X 2 )TtH 2 ] dH. 

For small values of X the distribution W involves two main types of contributions. The 
first arises when the change in the perturbation induces an anticrossing of the level v and 
can not be treated within perturbation theory. However, since it scales in proportion to the 
number of anticrossings encountered it vanishes as X^ +1 . The remaining contribution can 
be treated perturbatively in X. We will be concerned with the leading order of perturbation 
theory which becomes asymptotically exact in the limit X — > 0. The numerical simulations 
below show that even for small non-zero values of X this contribution gives an excellent 
approximation to W. To the first order of perturbation theory, an expansion around X 
gives (r — > a) 

dA„ 



where 



W(v, v; X) = (6(v - A V {X)) s(v - X-^ (X)) ), (5) 



dX 



E V {X) - E^X) 



and A U (X) = Q\ip v (a\ X)\ 2 denotes the intensity of the wavefunction at site a. Later we 
will discuss the role of the higher order contributions in controlling the asymptotic form of 
W and restricting the region over which Eq. (|[) is applicable. 

After Fourier transforming, and averaging over the Gaussian perturbation Hi, we obtain 



'" ivt 



W(v,v;x)=l —e- lvt W(v,t;X), (7a) 

-co Zll 



W(v, t; x) = ( 5(v — A v ) exp 



±2 2 AyA^ 

t x 



e lt - e v 



|2 



where we have made use of the rescaling with C(0) = X 2 /NA 2 



(7b) 
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Existing studies, such as the supersymmetry approach, show that, in the universal (long- 
time) limit, the statistical properties of the wavefunctions become decoupled from the spec- 
trum |LH]. As a result, the former can be averaged separately and we obtain Eq. (0), with 




00 dk 
-00 2n 



e~ isk B (k), (8a) 
(2/(3)(k/2nf]-^\ (gb) 



The remaining average in Eq. ( pb[) is performed over the Wigner-Dyson distribution of 
the N eigenvalues {e M }. 

The eigenvalue e u is identified by the correlator as special. In the limit iV — > 00 the 
translational invariance of the average allows us to set this eigenvalue to zero. This introduces 
an additional factor of | det H&\ into the average over the remaining N — 1 eigenvalues, and 
leads to the expression 

~ _ 1 / det H 2f3 \ 

p[ ) ~ \det[tf 2 + (2/(3)(kA/2rc)Y /2 / ' [ ) 

where Z N/3 = (\ detH^\) denotes a normalization constant. 

For unitary ensembles, an expression for Eq. (^|) can be found explicitly \Wj . The method, 



which is a generalization of an approach introduced by Guhr []17|| , relies on a super algebraic 
construction in which the spectral determinant is expressed as a Gaussian integral over 
fermion and bosonic variables. By integrating out the original degrees of freedom involving 
the matrix Hamiltonian H, an expression is found for Bp(k) which involves only a 6 x 6 
supermatrix. In contrast to the supersymmetry approach, instead of deriving a non-linear 
cx-model from the effective theory, we make use of a generalization of the Itzykson-Zuber 
integral [|l8l over a pseudounitary group to obtain the expression shown in Eq. @. 



Let us now consider the asymptotics of the distribution function. Eqs. (|) and (^) 
shows that large fluctuations arise when levels become anomalously close. In this limit 
v ~ (x\/jj/uj), where ui = \e u — e u ±\\ <C 1 denotes the smallest separation of two neighboring 
levels. RMT implies a probability distribution for level separations of p(u>) ~ uj 13 [l0[, which 

suggests 

W(v,v;x) ~ P p (v) p(oo)^ ~ Pp(v)^7= ^ +2 \ (10) 



with s = v/(2nx\/d). This agrees with the large s asymptotics in Eq. (3). As s increases the 
perturbative expansion eventually fails. In particular, for v > 1 higher orders of perturbation 
theory in x become important when s ~ while for v < 1 non-perturbative contributions 
begin to dominate at s ~ yti/x. At the end-points \v\ = 2v the distribution is essentially 
controlled by level anticrossings and Eq. (§) does not apply. One can show that Wpiv^v — > 
±2v- x) ~ xP +1 \v T 2v\ /3 / 2 - 1 P (3 (2v). 

The most striking feature of these results is the powerlaw decay of the correlator Bp(s), 
in contrast to the usual exponential behavior of other types of wavefunction distributions 
|16| . This admits to the existence of large fluctuations of the wavefunction even for small 



perturbations. The typical scale of the relaxation increases is proportion to x and y/v. 
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To verify the validity of the theoretical predictions, we have compared Eqs. (0) and ([3D 
with data taken from numerical simulations for the three symmetry classes. Each ensemble 
was defined by 45 realizations of a 7x 7x 7 tight-binding Anderson Hamiltonian. The diffusive 
behaviour of the majority of states was ensured by choosing the "on-site" disorder potential 
randomly from the range [—3,3], and a nearest-neighbor hopping integral normalized to 
unity. 

For the unitary case, we adopted the geometry of a torus for the lattice and an Aharonov- 
Bohm (A-B) flux was employed as the perturbation, X. To ensure unitary symmetry for 
even vanishing flux, fixed quasi-periodic boundary conditions were applied in the transverse 
directions. In both the orthogonal and symplectic cases, the boundary conditions in one 
direction were left open creating a lattice geometry of a ring. In both "radial" 
electric field served as the external parameter, X. Symplectic symmetry was imposed by 
spin-dependent hopping elements mimicing the action of a spin-orbit coupling. 

The numerical results were accumulated by taking an averaging over a wide range of 
values of the external perturbation, X. To minimize the systematic variation in A and C(0) 
as a function of energy, only a fifth of the total number states, located near the center of 
the band, were selected. 

The results for each Dyson ensemble are shown in Figs. |IJ ^|, and ^. In all cases, the 
factorization and scaling of Eq. (0) is born out. For v = 0, good quantitative agreement 
with the results implied by the \ 2 distribution is found for all three ensembles. Moreover, 
the results in Fig. |l|(a) give excellent agreement with the analytical expression in Eq. @. 
For the latter, the deviation from the theoretical prediction anticipated as a breakdown of 
perturbation theory, seems to lie beyond the statistical accuracy of this numerical work. Al- 
though we presently lack a complete analytical expression for -Bi(s) and B^s), the numerical 
data are consistent with the powerlaw decay predicted in Eq. (|TTJD . 

In this letter we have examined the sensitivity of wavefunction intensities in quantum 
chaotic systems to changes in an external perturbation. The rescaling implied by Eq. (0) has 
been verified numerically for all three symmetry classes and provides a concrete prediction 
for experiment. An analytical expression has been provided for unitary ensembles which 
shows a good quantitative agreement with numerical simulations. Our results indicate that 
variations of an external perturbation cause large wavefunction amplitude fluctuations to 
decay only as a power law. 

We wish to thank B. L. Altshuler, T. Guhr, C. H. Lewenkopf, F. von Oppen, and N. 
Taniguchi for valuable discussions. (BDS) and (AVA) wish to acknowledge the support of 
JSEP No. DAAL 03-89-0001, EC grant No. SCC-CT 90-00s20, and the hospitality of NECI 
where some of this work was performed. (BDS) is also greatful to the support of the Royal 
Society. 
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FIGURES 



FIG. 1. Comparison of the theoretical estimate of the distribution function W(v,v;x) with 
results taken from 45 realizations of a 7x 7x 7 Anderson model with an A-B flux (unitary ensemble), 
(a) shows the variation of W(v, v; x) with v for fixed values of x and v; (b) shows x x W(v, 0; x) 
for different values of x. The solid lines are the theoretical predictions. The inset shows the 
distribution of wave function intensities compared to the \ 2 distribution P2(s). 

FIG. 2. Same as in Fig. 1, but without A-B flux or spin-orbit coupling (orthogonal ensemble). 
The dashed lines in (a) and (b) are the result of fittings based on Eqs. (10) and (2), respectively. 
The inset shows the distribution of wavefunction intensities compared to Pi(s). 

FIG. 3. As in Fig. 2, but now including a spin-dependent hopping term (symplectic ensemble). 
The inset shows the distribution of wavefunction intensities compared to -Ri(s). 



7 




1 




2 



T3 

CD 
N 

• F— i 



'l-l 

o 

CD 
f-i 



10 l 



10" 



10" 



b~ 10 



10" 



-4 



(a) 







T 3H 
o 



-2 



\ o £ = 0.03 u = 4 
; □ £ = 0.08 u = 2 
* o i = 0.18 v = l 



5 = 



2 

17 / (2irXy/& 



CD 



X 



10 



-1 



10 



10 



-2 



10 



-4 



10 



-5 



10 u 



10" 



10"' 



10"' 



O £ = 0.03 
□ £ = 0.08 
O £ = 0.18 




10" 



10" 1 10° 



10 



o 



(b) 



3 



